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We present the results of Monte Carlo mass-loss predictions for massive stars covering a wide 
range of stellar parameters. We critically test our predictions against a range of observed mass- 
loss rates - in light of the recent discussions on wind clumping. We also present a model to 
compute the dumping-induced polarimctric variability of hot stars and we compare this with 
observations of Luminous Blue Variables, for which polarimetric variability is larger than for 
and Wolf-Rayet stars. Luminous Blue Variables comprise an ideal testbed for studies of wind 
clumping and wind geometry, as well as for wind strength calculations, and we propose they 
may be direct supernova progenitors. 



1 Introduction 



This contribution consists of two complementary as- 
pects of hot-star winds. We first describe the re- 
sults of mass-loss predictions - widely used in cur- 
rent massive star models in the galaxy, and beyond. 
In particular, we test our predictions as a function 
of effective temperature against recent radio data, 
and we discuss the potential implications for the 
clumping properties of supergiants of various spec- 
tral types (Sect. 2). We also discuss mass-loss pre- 
dictions for the winds of Luminous Blue Variables 
(LBVs), and we present results of the dumping- 
induced polarimctric variability of hot-star winds 
(Sect. 3), before we conclude. 



2 Monte Carlo mass-loss 
predictions 



2.1 Results: successes 

To gauge the success of our models and to be able 
to identify discrepancies, we test our predictions 
against a survey of radio mass-loss rates (Benaglia 
et al. 2007) from the free-free emission in hot-star 
winds. 
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Our method to predict the mass-loss rates of mas- 
sive stars is based on Monte Carlo radiative transfer 
calculations (sec Vink ct al. 1999, dc Koter et al. 
1997). In short, we compute non-LTE level pop- 
ulations for all relevant ions from hydrogen to zinc, 
before we follow the fate of a large number of photon 
packets through the wind. We predict the efficiency 
of the momentum transfer from the photons to the 
gas, generally assuming a pre-described velocity law 
(but see Vink et al. 1999, Muller & Vink, in prep.). 
We derive wind efficiencies, r\ = (Mu o0 )/(L/c), for a 
range of stellar parameters (including mctallicitics). 



Figure 1: Radio wind efficiency vs. effective tem- 
perature (Benaglia et al. 2007). Over- 
plotted are mass-loss predictions (Vink et 
al. 1999). Note the possible bi-stability 
jump at 21 000 K. 



In Fig. 1, we plot the wind efficiency versus ef- 
fective temperature in the regime of the wind bi- 
stability, where winds are predicted to change from 
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lower M, fast winds, on the hot side, to higher M, 
slow winds on the cool side - a result of a change 
in the Fe ionisation that drives the wind. We over- 
plot the mass-loss predictions around the bi-stability 
jump (dashed line) and focus on the general trends, 
before we continue our discussion on quantitative as- 
pects, and their implications, in Sect 2.2. 

The overall behaviour shows that r\ declines when 
the temperature drops. At the highest tempera- 
tures, the flux and the opacity show large overlap 
and the momentum transfer efficiency is maximal. 
At lower effective temperature, the flux is gradually 
emitted towards lower wavelengths, and there is a 
growing mismatch between the flux and the ultravi- 
olet line opacity. Figure f shows that our predictions 
of this overall behaviour are confirmed by the radio 
survey. Our second prediction is that there should 
be an increase in the mass-loss rate due to the opac- 
ity increase when Fe recombines from Fe iv to Fe 
hi. Our radio data appear to confirm the presence 
of a local maximum around 21 000 K, although data 
around this critical temperate is, as yet, sparse. 

The bi-stability limit is relevant for stellar evo- 
lution calculations when stars evolve off the main 
sequence towards the red part of the Hertzsprung- 
Russell diagram (HRD). It may also play a role for 
LBV winds, when LBVs such as AG Car, change 
their temperatures - and radii - on timescales of the 

5 Dor variations (of the order of years to decades). 
This variable wind behaviour - predicted by Vink 

6 de Koter (2002) - is anticipated to result in cir- 
cumstcllar media consisting of concentric shells with 
varying wind densities. Kotak & Vink (2006) re- 
cently suggested that the quasi-periodic modulations 
seen in the radio lightcurves of some supernovae may 
imply that LBVs could be direct supernova progen- 
itors. At first this seems at odds with stellar evo- 
lution calculations, which do not predict LBVs to 
explode. However, there is a growing body of ev- 
idence hinting that LBVs may nonetheless explode 
(Pastorello et al. 2007, Smith 2007, Gal- Yam et al. 
2007). 

Despite the success of our models in explaining 
LBV mass-loss variability, the bi-stability jump, and 
the scaling of M with metallicity (see de Koter, 
this volume), we turn to discrepancies of our models 
against empirical mass-loss rates. 

2.2 Results: discrepancies 

Discrepancies have been noted between the Vink et 
al. (2000) predictions and empirical mass-loss rates 
in several areas of the HRD. One group of objects is 
that of low L (log L/Lq < 5.2) O dwarfs, where the 
data of Martins et al. (2005) fall well below predic- 
tions, by factors of 10 or more. The reason for this 
discrepancy is as yet not understood. Another area 
of discrepancies is that of the B supcrgiants where 
empirical rates have been found to be much lower 



than predicted rates (Vink et al. 2000, Trundle & 
Lennon 2005). 

The most worrisome however is the situation with 
garden- variety type O stars! Figure 1 shows that the 
Vink et al. (2000) predictions are lower than the ob- 
served rates. The radio rates are likely to be upper 
limits as the radio free-free emission is a p 2 process, 
and any form of clumping leads to a maximal M. 
Mokiem et al. (2007) and Puis (this volume) noted 
that if this discrepancy is related to wind clumping 
- and theoretical rates are unaffected - the empirical 
p 2 mass-loss rates must be down-revised by a factor 
2-3, suggesting a clumping factor / ~5. Recent mas- 
sive star evolution models would not be effected by 
such modest clumping factors as these already use 
the theoretical Vink et al. (2000) rates. 

If / ^5 were universal, one would expect the em- 
pirical M for B supergiants to be lower than indi- 
cated in Fig. 1 and the discrepancy would amount to 
an order of magnitude, or more. This implies there 
are some serious issues with our theoretical under- 
standing of hot-star winds, and we need to recon- 
sider even our most basic modelling assumptions, 
such as sphericity and homogeneity, which can be 
tested with linear polarimctry. 



3 Linear polarisation variability 
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Figure 2: Polarisation variability vs. stellar radius 
- for different clump ejection rates. N = 
N ifl, where N is the clump ejection rate 
(related to the mass-loss rate) and t$. = 
R/voo (see Davies et al. 2007). 



Linear polarimetry is a tool to measure asymmetries. 
Davies et al. (2005) performed a polarimetric sur- 
vey of LBVs and found asymmetries in a majority 
of them. When the position angle (PA) of the polar- 
isation shows a straight line in the Stokes QU dia- 
gram this is generally attributed to a large-scale, axi- 
symmctry, e.g. a disk. Davies et al. (2005) found 
time-variable PAs for objects such as AG Car and 
attributed these to wind clumping. Subsequently, 



2 



Advances in mass-loss predictions 



Davies ct al. (2007) constructed an analytic clump- 
ing model, releasing clumps with a certain ejection 
rate per wind flow-time, TV, from the wind base. The 
average polarisation of the clump ensemble was cal- 
culated; the results for LBVs with a range of tem- 
peratures and radii are shown in Fig. 2. 

When the LBVs decrease their radii, the clumps 
become smaller and denser, and produce more po- 
larisation. This behaviour reverses at the tempera- 
ture of the bi-stability jump where the wind becomes 
faster, and the clumps spend less time at the wind 
base. As a result the polarisation drops. Figure 2 
also shows that when the ejection rate increases, the 
polarisation drops as the wind approaches that of a 
smooth outflow leading to zero polarisation. 



~ 2 





i >- 


1 1 1 1 1 1 " 

* 







V 


« 

_ • 




\ 

• 

V 

• 


P Cyg observed . 
i i i i i . . . 


0.1 


1.0 


10.0 100.0 1000.0 



w 

Figure 3: Polarisation variability vs. clump ejection 
rate. The presence of a maximum implies 
there are two solutions compared to the 
observed level of P Cyg. 

We now consider the average polarisation as a 
function of ejection rate. The right-hand side of 
Fig. 3 shows the regime of many optically thin 
clumps. A maximum is reached at TV ~10 where 
the clumps become optically thick and multiple scat- 
tering becomes important. The left-hand side rep- 
resents the optically-thick clump regime. The ob- 
served level of polarisation of P Cyg is shown as a 
horizontal line that intercepts both branches. The 
data either point to a wind with low TV, or to one 
with TV ~ 1000. We distinguish between these two 
branches using timescale information of the polari- 
metric variations. Preliminary results from our re- 
cent monitoring campaign indicate that the high- TV 



scenario is likely to be the correct one (Davies ct al. 
in prep) , which would suggest that LBV winds con- 
sist of thousands of optically thin clumps close to 
the photosphere. 

4 Conclusion 

Polarimetry is a tool to constrain the clumping prop- 
erties of hot-star winds. This may become a pow- 
erful means by which to constrain non-LTE mod- 
els and mass-loss predictions. We anticipate to wit- 
ness an increased understanding of hot-star winds - 
an important endeavour because of the impact mass 
loss has on massive star evolution modelling. 
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